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Abstract 

We explore the effects of nuclear masses on the temperature and neutron density conditions 
required for r-process nucleosynthesis using four nuclear mass models augmented by the latest 
atomic mass evaluation. For each model we derive the conditions for producing the observed 
abundance peaks at mass numbers A ~ 80, 130, and 195 under the waiting-point approximation 
and further determine the sets of conditions that can best reproduce the r-process abundance 
patterns (r-patterns) inferred for the solar system and observed in metal-poor stars of the Milky 
Way halo. We find that (1) the conditions for producing abundance peaks at A ~ 80 and 195 
tend to be very different, which suggests that these two peaks are not produced simultaneously; 
(2) the typical conditions required by the critical waiting-point (CWP) nuclei with the N = 126 
closed neutron shell overlap significantly with those required by the N = 82 CWP nuclei, which 
enables coproduction of abundance peaks at A ~ 130 and 195 in accordance with observations of 
many metal-poor stars; and (3) the typical conditions required by the N = 82 CWP nuclei can 
reproduce the r-pattern observed in the metal-poor star HD 122563, which differs greatly from the 
solar r-pattern. We also examine how nuclear mass uncertainties affect the conditions required 
for the r-process and identify some key nuclei including 76 Ni to 78 Ni, 82 Zn, 131 Cd, and 132 Cd for 
precise mass measurements at rare-isotope beam facilities. 
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I. INTRODUCTION 



Nucleosynthesis via rapid neutron capture, the r-process, is a major mechanism for pro- 
ducing the elements heavier than Fe [HQ- Understanding this process requires knowledge 
of properties such as masses, /3-decay lifetimes, and neutron-capture cross sections for a 
large number of extremely neutron-rich nuclei far from stability (e.g., [35]). Most of this 
nuclear input is beyond the reach of experiments in the foreseeable future and therefore, 
must be calculated with guidance from existing data and from measurements to be made at 
rare-isotope beam facilities such as FRIB, RIKEN, and FAIR. In this paper we explore the 
importance of nuclear masses in estimating the temperature and neutron density conditions 
required for a specific r-process scenario, where neutron-capture reactions are in equilibrium 
with the reverse photo-disintegration reactions, i.e., there is (^,7) ^ (7,^) equilibrium 
(e.g., |3rl5|). Using four nuclear mass models, we show that the required conditions can be 
determined mostly from the neutron-separation energies for a small number of critical nuclei 
with iV = 50, 82, and 126 closed neutron shells. For each model, we further determine the 
best-fit sets of conditions to reproduce the r-process abundance pattern (r-pattern) inferred 
for the solar system and observed in metal-poor stars of the Milky Way halo. This allows 
us to draw several interesting conclusions regarding the production of different parts of the 
overall r-pattern in an (n, 7) # (7, ^-equilibrium scenario. We also illustrate the effects 
of nuclear mass uncertainties on the required r-process conditions and identify the key nu- 
clei that have the largest impact and therefore, are important candidates for precise mass 
measurements at rare-isotope beam facilities. 

We first give a brief overview of the r-process. Detailed reviews can be found in Ref. {3^ 
g Historically, the abundance distribution of nuclei in the sola, systen, played a ccncM 
role in studies on the origin of the elements 12J. One of the prominent features of this 
distribution is the presence of three sets of double peaks in the region beyond the Fe group 
nuclei. This was recognized as signatures of two distinct processes of neutron capture: a slow 
(s) one encountering the N = 50, 82, and 126 closed neutron shells in the stable region and 
a rapid (r) one encountering the same in the extremely neutron-rich region of the nuclear 
chart [l|, Q]. Specifically, the peaks at mass numbers A ~ 80, 130, and 195 were produced 
by the r-process and represent the crucial features of the solar r-pattern, which is derived 
by subtracting the s-process contributions from the net solar abundances (e.g., [6j). 
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In order to fully understand the r-process, we need conditions such as temperature and 
neutron density in the associated astrophysical environments in addition to the properties 
of a large number of extremely neutron-rich nuclei. Neither the astrophysical nor the nu- 
clear input is firmly established, although much progress has been made over the past two 
decades {3- 5]. Proposed astrophysical sites for the r-process include neutrino-driven winds 
from proto-neutron stars formed in core-collapse supernovae (CCSNe) [z-9], shocked surface 



layers of O-Ne-Mg cores associated with low-mass CCSNe 
of black holes formed in high-mass CCSNe 



and ejecta from neutron star mergers 
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101 ]. winds from accretion disks 



14j . He shells of metal-poor CCSNe 
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20] . There are large uncertainties in the conditions 



associated with all CCSNe environments due to the substantial uncertainties in modeling 



such environments (e.g., [21 
matter is considered (e.g. 



22 



especially when neutrino transport in hot and dense nuclear 



23|). Whi 



e recent studies lend much support to neutron 



star mergers being an r-process site 17H20I]. it remains to be seen whether such models are 



consistent with the history of r-process enrichments in the Milky Way and in its satellite 



dwarf galaxies (e.g. 



24 



Further, how sensitive these models are to the uncertainties in 



the current understanding of the nuclear equation of state (e.g., (27]) and to the numerical 
treatment of the merger dynamics remains to be studied in detail. 

The conditions in the astrophysical environments relevant for the r-process ultimately 
boil down to the seed nuclei for neutron capture at the beginning of the process and the 
temperature T(t ) and the neutron (number) density n n (t) as functions of time t during the 



process (e.g. 



28||). In cases where neutrino interactions are important, the time evolution of 



neutrino fluxes and energy spectra is also required (e.g., 15|, ll6|]). In the rest of the paper we 
will ignore neutrinos and focus on a broad class of astrophysical environments where matter 
undergoes r-processing at T > 10 9 K and n n > 10 20 cm -3 . For such high temperatures and 
neutron densities, previous studies have shown that (n, 7) (7, n) equilibrium is achieved 



[e.g., 



29]). In this equilibrium, the abundance distribution in each isotopic chain at a specific 



proton number Z is almost strongly peaked at one nucleus. This is referred to as a waiting- 
point (WP) nucleus because upon reaching it, the r-process must wait for it to /3-decay 
before producing heavier nuclei. Under this so-called WP approximation, the r-process 
path is defined by all the WP nuclei heavier than the seed nuclei, and the progress along 
this path is regulated by the /3-decay of these WP nuclei. So long as this approximation is 
valid, there is no need to follow neutron capture and photo-disintegration reactions, which 
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greatly simplifies the r-process calculation. 



Due to the equilibrium between neutron capture and photo-disintegration reactions, the 
abundance ratio between two neighboring isotopes is given by the Saha equation (e.g., [3|-|5j]): 



Y(Z,A + 1) 



2nh 



2X§ 



n. 



G(Z,A+l) /A + iy 



exp 



S n (Z,A+l) 
kT 



1) 



Y(Z,A) \m u kT J 2G(Z,A) V A 

where % is the Planck constant, m u is the atomic mass unit, k is the Boltzmann constant, 
(Z, A) indicates a nucleus with proton number Z and mass number A, and Y, G, and 
S n denote the number abundance, partition function, and neutron separation energy of 
the appropriate nucleus, respectively. For a specific isotopic chain, the corresponding WP 
nucleus has the largest abundance and is determined by the partition functions and neutron 
separation energies of the relevant nuclei for fixed T and n n . As can be seen from the 
exponential dependence on the neutron separation energy in Eq. (JT]), nuclear masses are 
among the most important input for modeling the r-process. The other crucial input is 
/3-decay lifetimes of the relevant nuclei. 

Over the past two decades, tremendous progress has been made in measuring nuclear 
properties relevant for the r-process. For example, the /3-decay half-lives of 38 very neutron- 
rich isotopes bordering the r-process path have been measured recently j^]. In addition, 
the masses of a group of nuclei including 80 Zn 3l|, |32j and 130 Cd 33J have been measured 
with a very high accuracy 34|. Meanwhile, considerable advance has been made in the 
theoretical investigation of nuclear masses. The four nuclear mass models used in this paper 
span from the macroscopic-microscopic kind, represented by the finite-range droplet model 
(FRDM) 35| and a more recent Weizsacker-Skyrme (WS*) model 36j, to the microscopic 
kind, represented by the Skyrme-Hartree-Fock-Bogolyubov mean-field (HFB-17) model 37] 
and the relativistic mean- field (RMF) model [38?]. These models can reproduce the exper- 
imentally known neutron separation energies with a root-mean-square deviation of 0.399 
(FRDM), 0.332 (WS*), 0.506 (HFB-17), and 0.653 (RMF) MeV, respectively. 

Based on the above overview, there are two frontiers of r-process research: one focusing 
on the search for the astrophysical sites and quantification of the conditions therein and 
the other on acquiring a reliable database for the relevant nuclear input. Observations of 
elemental abundances in metal-poor stars of the Milky Way halo (see [39| for a review) have 
shed important light on the r-process sites (e.g., 40|, |41|)- The observed r-patterns also 
provide an important test of the basic soundness of the nuclear input (e.g., [42|). Of course, 



the astrophysical and the nuclear input must be coupled together in order to produce an 
r-pattern for comparison with observations. With substantial uncertainties in the current 
understanding of both the r- process sites and the nuclear input, parametrization of the 
astrophysical conditions is often used in exploring the effects of nuclear input on the r-process 



production (e.g., 



28 



421]). As a practical matter, we adopt the classical approach of using 



T, n n , and the corresponding neutron irradiation time r along with the WP approximation 



(e.g., 



431 ] ) to carry out our r-process calculations below. Our main purpose is to explore 



the effects of the four nuclear mass models mentioned above on the T and n n conditions 
required for r-process nucleosynthesis. 

We give a detailed discussion of the WP approximation in Sec. [TTJ Using this approxima- 
tion along with four nuclear mass models, we derive for each model the T and n n conditions 
that are required for producing the abundance peaks at A ~ 80, 130, and 195 as observed 
in the solar system. In Sec. IIHI we describe the classical approach to simulate the r-process 
and use this approach to determine the sets of conditions that can best reproduce the solar 
r-pattern and the r-patterns observed in metal-poor stars for each of the adopted nuclear 
mass models. We discuss our results and give conclusions in Sec. [TV] 



II. r-PROCESS CONDITIONS UNDER THE WP APPROXIMATION 

While there are substantial uncertainties in both the r-process sites and the relevant 
nuclear input, an essential feature of the r-process is considered robust: the observed abun- 
dance peaks at A ~ 80, 130, and 195 correspond to the intrinsic properties of extremely 
neutron-rich nuclei with N = 50, 82, and 126 closed neutron shells that are produced in 
the r-process. As discussed in Sec. HI when (n, 7) ^ (7, n) equilibrium is achieved, the 
total abundance of an isotopic chain is concentrated in the corresponding WP nucleus. The 
/9-decay lifetimes of the WP nuclei then regulate the abundance pattern resulting from an 
r-process episode. In particular, the much longer /3-decay lifetimes of extremely neutron-rich 
nuclei with closed neutron shells than those without produce peaks in r-patterns. Conse- 
quently, in order to produce the observed peaks in r-patterns under the WP approximation, 
nuclei with N = 50, 82, and 126 closed neutron shells and with A ~ 80, 130, and 195, re- 
spectively, must be among the WP nuclei. These critical WP (CWP) nuclei were identified 
by earlier studies (e.g., [44j). Table [I] lists the CWP nuclei in the present work. These nuclei 
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are selected by referring also to previous r-process calculations 44J. Below we follow the 
spirit of these studies to derive the conditions required for r-process nucleosynthesis under 
the WP approximation by considering the properties of the CWP nuclei. 



TABLE I: Critical waiting-point (CWP) nuclei 

N CWP nuclei 

50 80 Zn, 79 Cu, 78 Ni 

82 130 Cd, 129 Ag, 128 Pd, 127 Rh, 126 Ru 

126 195 Tm, 194 Er, 193 Ho, 192 Dy, 191 Tb 



A. Role of neutron separation energies 

To be quantitative, we define a WP nucleus (Z, A^vp) as one that has an abundance 

Y(Z,A WP ) >0.5J2Y(Z,A), (2) 



.4 



where the sum over A gives the total abundance of the corresponding isotopic chain. For 
specific T and n n , we can use the above criterion and the relative abundance Y(Z, A + 
1)/Y(Z, A) given by Eq. (OQ) to determine (Z, Awp) from the nuclear partition functions and 
neutron separation energies provided by a model. Conversely, we can also determine the 
T and n n conditions required by a specific WP nucleus. As can be seen from Eq. (CQ), the 
predominant dependence of Y(Z, A + 1)/Y(Z, A) is on the neutron separation energy, which 
can be calculated from a nuclear mass model. We ignore the small differences in the nuclear 
partition function and in the mass number and rewrite Eq. (JT|) as 



Y(Z,A + 1) 
Y(Z,A) 



exp 



S n (Z,A + l)-SZ(T,nr, 
kT 



(3) 



where 
S° n (T,n n ) = kT In 



2 



m u kT^ 
2tt?i 2 I 









^ 3/2- 


= T 9 


2.79 + 



1.5 log T 9 -log(n n /10 20 cm 



-3\ 



5.04 



MeV. 



(4) 



In the second equality of Eq. (j3J), T 9 is T in units of 10 9 K. Equation ((3]) is used in the 
calculations below. 

As an example, we adopt the WS* mass model to calculate the sets of Tg and n n within 
the ranges 1 < T 9 < 3 and 10 20 < n n < 10 30 cm" 3 that are required by the N = 82 CWP 
nuclei. The results are shown in Fig. [TJ For a specific Tg, the values of n n between two 
identical symbols in this figure would allow the corresponding nucleus to have > 50% of the 
total abundance of its isotopic chain. In order to accommodate all the iV = 82 CWP nuclei, 
the common range of n n for a specific Tg is bounded from below by 126 Ru (filled circle) and 
from above by 130 Cd (filled triangle). This range of n n changes with T 9 and is shown as the 
shaded band in Fig. [TJ This band represents the T 9 -n n conditions required by the N = 82 
CWP nuclei. 
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FIG. 1: The Tg-n n conditions required by the N = 82 CWP nuclei based on the WS* mass model. 
For a specific Tg, the values of n n between two identical symbols would allow the corresponding 
nucleus to have > 50% of the total abundance of its isotopic chain. The conditions indicated by 
the shaded band are required to accommodate all the N = 82 CWP nuclei. 



While the neutron separation energies of a large number of nuclei from the WS* mass 
model are used in the above calculations, the results in Fig. [TJ are determined effectively by 
the two-neutron separation energies of the N = 82 CWP nuclei and their N = 84 isotopes 
due to nuclear systematics. This can be understood as follows. Because of the effect of 
pairing on neutron binding, all WP nuclei have even N. The relative abundance of two 
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neighboring even-iV isotopes can be obtained from Eq. ([3]) as 

~S 2n (Z,A)-2S°(T,n r . 



Y(Z,A) _ 
Y(Z,A-2)~ eXP 



kT 



(5) 



where S 2n denotes the two-neutron separation energy. The values of S 2n /2 for isotopes of 
Ru, Rh, Pd, Ag, and Cd (44 < Z < 48) are shown as functions of N in Fig. 121 which exhibits 
the general trend that S 2n essentially monotonically decreases with iV for a specific isotopic 
chain. Based on this aspect of nuclear systematics, it can be seen from Eq. ([5]) that the abun- 
dance of an even-iV isotope increases with N [Y(Z, A)/Y(Z, A — 2) > 1] until S 2n (Z, A)/2 
falls below S®(T,n n ), from which point on it decreases with N [Y(Z, A)/Y(Z, A — 2) < 1]. 
Therefore, the abundance of even-iV isotopes peaks at the nucleus (Z,Awp), for which 

S 2n (Z,A WP + 2) < S° n (T,n n ) < S 2n (Z,A WP ). (6) 

The above equation effectively defines a WP nucleus (e.g., j45() and can be used to determine 
the Tg-n n conditions required by a specific WP nucleus. For example, 130 Cd has S 2n /2 = 



5.488 MeV while 132 Cd has S 2n /2 = 2.869 MeV [36j. So the T 9 -n n conditions corresponding 
to 2.869 < S°{T,n n ) < 5.488 MeV [the band between dashed lines labeled as S°(Cd) in 
Fig. |2] are required for 130 Cd to be a WP nucleus. Likewise, the shaded band labeled as 
S® in Fig. [2] corresponds to the conditions required to accommodate all the N = 82 CWP 
nuclei (shaded band in Fig. [p. 



B. Tg-n n conditions for four nuclear mass models 

The calculations in Sec. Ill Al can be generalized to determine the Tg-n n conditions required 
by the iV = 50, 82, and 126 CWP nuclei, respectively, for any specific nuclear mass model. 
The results for the N = 50 CWP nuclei are presented for the FRDM, WS*, and RMF models 
in Fig. |3k, which clearly show that the required conditions change with models. Similar to 
the case of the iV = 82 CWP nuclei (Figs. [1] and |2]) discussed in Sec. Ill the upper curve 
for each model in Fig. |3^ is effectively determined by the two-neutron separation energy of 
the lightest iV = 50 CWP nucleus 78 Ni and the lower curve by that of the N = 52 isotope 
82 Zn of the heaviest N = 50 CWP nucleus 80 Zn. Therefore, the large differences among the 
conditions required by the iV = 50 CWP nuclei for different models can be traced to the 
differences in the two-neutron separation energies of 78 Ni and 82 Zn provided by these models. 
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FIG. 2: (Color online) Two-neutron separation energies S% n for even-iV isotopes of Ru, Rh, Pd, 
Ag, and Cd (44 < Z < 48) in the WS* model shown in terms of Si n j2 as a function of iV around 
N = 82. The band between the dashed lines labeled as 5°(Cd) corresponds to the Tg-n n conditions 
required for 130 Cd to be a WP nucleus. Likewise, the shaded band labeled as S® corresponds to 
the conditions required to accommodate all the N = 82 CWP nuclei (shaded band in Fig. [I]). See 
text for details. 

In particular, the differences for 82 Zn among the models appear to be substantially larger 
than those for 78 Ni. We also note that no conditions can be found to accommodate all the 
N = 50 CWP nuclei for the HFB-17 model, for which the odd-even effects in the neutron 
separation energy for Ni, Cu, and Zn isotopes around N = 50 are larger by ~ 1-1.5 MeV 
than those for the FRDM, WS*, and RMF models. 

Noting that the two- neutron separation energies of 78 Ni and 82 Zn can be calculated from 
the masses tabulated in the latest atomic mass evaluation AME2011-preview 3^], we aug- 
ment the nuclear mass models by using the tabulated values in AME2011-preview when they 
are available to replace the corresponding model predictions. Remarkably, all four models, 
including the HFB-17 model, now give the same conditions required by the A^ = 50 CWP 
nuclei as shown in Fig. [3b- We find that the changes between Figs. [3k and [3b are caused 
dominantly by the use of the tabulated masses of 76 Ni to 78 Ni and 78 Zn to 82 Zn, which con- 
firms the crucial roles of the two- neutron separation energies of 78 Ni and 82 Zn in determining 
the conditions required by the A^ = 50 CWP nuclei. In the calculations below, we use the 
FRDM, WS*, HBF-17, and RMF models that are augmented by AME2011-preview. 

We calculate the conditions required by the A^ = 82 and 126 CWP nuclei, respectively, 
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FIG. 3: (Color online) The Tg-n n conditions required by the N = 50 CWP nuclei, (a) The band 
between two curves of the same kind represents the required conditions based on the corresponding 
mass model. Note the large differences among the results for the three indicated models. Note 
also that no conditions can be found to accommodate all the N = 50 CWP nuclei based on the 
HFB-17 model, (b) Same as (a), but for those nuclei the masses of which are tabulated in the 
latest atomic mass evaluation AME2011-preview [341 ] . the model predictions are replaced by the 
tabulated values. All four models now give the same conditions required by the A = 50 CWP 
nuclei. 

as in the case of the N = 50 CWP nuclei. The results are summarized in Fig. HI It can be 
seen that the conditions required by the A^ = 82 CWP nuclei (shaded band) are essentially 
converged for the four augmented nuclear mass models just like those required by the A^ = 50 
CWP nuclei (horizontally hatched band between solid curves). In contrast, the conditions 
required by the A" = 126 CWP nuclei (vertically hatched band between dashed curves) 
are still strongly dependent on models. This is because these nuclei and the majority of 
those in the nearby region of the nuclear chart are still out of the reach of experiments 
while theoretical predictions for their masses involve dramatic extrapolations with large 
uncertainties (e.g., j^). 




Figure H] resembles a phase diagram in terms of three bands for the production of the 
peaks at A ~ 80, 130, and 195 in r-patterns that correspond to the A" = 50, 82, and 126 
CWP nuclei. For the T$-n n conditions inside the non-overlap region of a band, only a single 
peak can be produced. For those conditions inside the overlap region of two bands, it is 
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FIG. 4: (Color online) The Tg-n n conditions required by the N = 50 (horizontally hatched band 
between solid curves), 82 (shaded band), and 126 (vertically hatched band between dashed curves) 
CWP nuclei, respectively, for four nuclear mass models: (a) RMF, (b) HFB-17, (c) FRDM, and 
(d) WS*. See text for details. 

possible to produce two peaks simultaneously. For the RMF model only, there is a very 
thin sliver where three bands overlap. Consequently, we consider it very unlikely that three 
peaks can be produced simultaneously. A close examination of Fig. H] shows that the Tg-n n 
conditions required by the N = 50 CWP nuclei (horizontally hatched band between solid 
curves) are distinct from those required by the N = 126 CWP nuclei (vertically hatched 
band between dashed curves) for the FRDM, HFB-17, and WS* models. These two sets 
of conditions overlap only slightly for the RMF model. This suggests that the peaks at 
A ~ 80 and 195 in r-patterns are not produced simultaneously. Their production may differ 
in the time of occurrence within the same astrophysical site or in the astrophysical site itself. 
In contrast, there is large overlap between the conditions required by the N = 50 and 82 
CWP nuclei (shaded band) for the four models considered. In addition, there is slight to 
significant overlap between the conditions required by the N = 82 and 126 CWP nuclei for 
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all the models. Therefore, it is possible to produce the peaks at A ~ 80 and 130 or those at 
A ~ 130 and 195 simultaneously. The above results will be examined by detailed r-process 
calculations in Sec. IIHI 



III. r-PATTERNS FROM THE CLASSICAL APPROACH 

Our main goal here is to explore the effects of nuclear masses on the conditions required 
for r-process nucleosynthesis under the WP approximation. As discussed in Sec. HI1 these 
conditions are mostly set by the neutron separation energies, which we calculate from four 
nuclear mass models augmented by the latest atomic mass evaluation AME2011-preview. 
As confirmation of these results, we calculate the r-patterns produced under the conditions 
shown in Fig. H] and compare them with those inferred for the solar system and observed 
in metal-poor stars. As the range of conditions shown in Fig. H] is rather broad, we sample 
these conditions at a fixed temperature. Specifically, we use T 9 = 1.5 and n n = 10 2ao -10 22 ' 5 , 
10 20 ' 5 -10 25 , and 10 23 ' 5 -10 27 ' 5 cm" 3 as typical conditions required by the N = 50, 82, and 
126 CWP nuclei, respectively. As noted in Sec. Ill there is overlap between these sets of 
conditions, which can lead to coproduction of two peaks in the r-pattern. 



Sis 



as follows. 



We carry out an r-process calculation using the classical approach (e.g., [3Hf 
We take the seed nucleus to be 56 Fe (Z = 26). We assume that some material with an initial 
abundance Yp e of 56 Fe is irradiated with neutrons at fixed T and n n for a time r. We then 
use Eq. ([3]) along with a nuclear mass model to calculate the fractional abundance P(Z, A) 
for all the nuclei in each of the isotopic chains with Z > 26, where 

Note that for fixed T and n n , P(Z,A) is also fixed. Using the fractional abundances, we 
calculate the effective /3-decay rate of an isotopic chain as 

\^ Z = ^2P(Z,A)\^Z,A), (8) 

A 

where A / g(Z, A) is the /3-decay rate of the nucleus (Z,A). For all our calculations, we use 
\p{Z, A) from the experimental data in Ref. (46] for the nuclei with measurements and from 
the theoretical estimates in Ref. j^] based on the FRDM+QRPA method for those without. 
We then solve the set of equations 

Y z (t) = -\p,zYz(t), Z = 2Q, (9) 
13 



Y z (t) = \p,z-iY Z -i(t) - X M Y z {t), Z > 26, 



(10) 



where the dot denotes the derivative with respect to time t, and Y z (t) is the total abundance 
of the isotopic chain with proton number Z at time t. The initial conditions are Y z (0) = Yp e 
for Z = 26 and for Z > 26. We assume that the r-process freezes out instantaneously at 
t = t. The freeze-out abundance of the nucleus (Z, A) is 

Y {o (Z,A) = P(Z,A)Y z (r). (11) 

The final abundance distribution from an r-process episode is obtained by following the (3 
and a decays of all the nuclei in the freeze-out distribution. We include /3-delayed emission 
of up to three neutrons [47| ,which has the important effect of smoothing the final r-pattern. 

Fission is 



The data on a-decays are taken from the National Nuclear Data Center [48 ] 
ignored in all the calculations. Note that the WP approximation is implicit in the classical 
approach as Y z (t) is dominated by the corresponding WP nucleus with P(Z, A WP ) > 0.5. 

For comparison with the r-patterns inferred for the solar system and observed in metal- 
poor stars, we need to superpose the results from many r-process episodes described above. 
As we take T 9 = 1.5 for all the calculations, we denote each episode by its n n . The neutron 
irradiation time r(n„) and the weight co(n n ) for each episode are taken to be 

r(n n ) = axn b n , (12) 
u){n n ) = cxn d n1 (13) 

where a, b, c, and d are parameters to be determined by a least-squares fit to the r-pattern 
used for comparison. The above procedure has been used extensively in r-process studies 



(e.g., 



28 



43 



49|-|52|). 



A. Comparison with solar-like r-patterns 

We first use the classical approach to reproduce the solar r-pattern for 125 < A < 209 



53[ (see also e.g., 5jJ, |54j), which is shown as the dashed curve in Fig. [5j We consider a 
superposition of nine neutron densities (equidistant on log scale) within the range 

10 23.5 < 

n n < 10 27 5 cm" 3 , which corresponds to the typical conditions required by the N = 126 
CWP nuclei for T g = 1.5 (see Fig. H]). The best-fit results (hereafter "Fit I") for the four 
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adopted nuclear mass models are shown as the solid curves in Fig. |5j Note that although 
the fits are performed for the solar isotopic r-pattern, the patterns shown in Fig. |5] are for 
the corresponding elemental abundances. It can be seen that the solar r-pattern from the 
peak at A ~ 130 (Z ~ 52, Te) to that at A ~ 195 (Z ~ 78, Pt) is reproduced rather well for 
the FRDM, HBF-17, and WS* models. For the RMF model, the rare-earth elements with 
Z = 66-70 (Dy, Ho, Er, Tm, Yb) are severely underproduced. This deficiency may reflect the 
necessity to adopt improved RMF parameter sets (e.g., PC-PK1 55j which provides much 
better descriptions for the nuclear ground-state and excited state properties [56j) or that 
the classical approach is inadequate to give a full description of r-process nucleosynthesis 
(e.g., instantaneous freeze-out is not a good approximation |28|). In any case, the peaks 
at A ~ 130 and 195 are reproduced adequately for all four mass models, which suggests 
that the typical conditions required by the N = 126 CWP nuclei can indeed produce both 
these peaks. As discussed in Sec. Ill B\ this is because there is significant overlap between 
the conditions required by the N = 82 and 126 CWP nuclei (see Fig. H|). 

Observations show that Ba [Z = 56) and heavier elements in many metal-poor stars of 
the Milky Way halo follow the solar r-pattern rather closely [39]. The values of loge(E 
log(E/H) + 12 for the elements with 38 < Z < 79 obtained by averaging the data 



54 



57 



39, 



58) on two such "r-II" stars, CS 22892-052 and CS 31082-001, are shown as filled 



circles in Fig. El The data on Pb (Z = 82) and Th (Z = 90) for CS 22892-052 are also 
shown. The solar r-pattern has been translated to pass through the data on Eu (Z = 63) 
and can be seen to represent the pattern for Z > 56 in r-II stars very well. Recently, Te 
(Z = 52) has been measured in a group of metal-poor stars (BD +17°3248, HD 108317, 



HD 128279 1591; HD 160617 



601]) . This extends the comparison of r-patterns in metal-poor 



stars with the solar r-pattern to include an element in the peak at A ~ 130. The data on 
Te for BD +17°3248 (shifted according to the observed Te/Eu ratio) is shown as the filled 
square in Fig. [5l It can be seen that the Te data is consistent with the solar r-pattern 
and with the coproduction of the peaks at A ~ 130 and 195 under the typical conditions 
required by the N = 126 CWP nuclei. The above results are also in agreement with previous 



studies (e.g., 49_|, |5l|), which concluded that the r-process responsible for the elements with 
56 < Z < 82 is characterized by neutron densities of 10 23 -10 28 cm -3 . 

The Fit I results shown in Fig. [5] cannot adequately reproduce the abundances of the 
elements with 38 < Z < 47 in r-II stars (especially when the RMF model is used). Further, 
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FIG. 5: (Color online) Elemental abundances (solid curves) calculated from the classical r-process 
approach using four nuclear mass models: (a) RMF, (b) HFB-17, (c) FRDM, and (d) WS*. The 
filled circles represent the average data on r-II stars of the Milky Way halo and the dashed curve is 
the solar r-pattern translated to pass through the Eu data. The filled square gives the Te abundance 
recently measured in the metal-poor star BD +17°3248 (shifted according to the observed Te/Eu 
ratio). The calculated solid curves (Fit I) are the best-fit results to reproduce the solar isotopic 
r-pattern for 125 < A < 209 with T Q = 1.5 and 10 23 - 5 < n n < 10 27 5 cm , which are the typical 
conditions required by the N = 126 CWP nuclei. See text for details. 

the elements in the peak at A ~ 80 (Z ~ 34) of the solar r-pattern are severely underpro- 
duced by these calculations. Additional r-process contributions to or alternative sources for 
the elements below the peak at A ~ 130 are thus required and this issue has been under 



active investigation 



40 



61 



64j. Here we explore the possibility that there are additional 
contributions from r-process nucleosynthesis under the conditions required by the N = 50 
CWP nuclei. We consider a superposition of six neutron densities (equidistant on log scale) 
within the range 10 2ao <n n < 10 22 ' 5 cm -3 to best reproduce the solar isotopic r-pattern for 
69 < A < 124 (hereafter "Fit II"). The results are shown as the solid curves in Fig. |6j It 
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can be seen that the conditions required by the N = 50 CWP nuclei indeed can produce the 
peak at A ~ 80. However, it is also clear that the r-patterns from the peak at A ~ 130 to 
that at A ~ 195 inferred for the solar system and observed in r-II stars require very different 
conditions from those for producing the peak at A ~ 80. 
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FIG. 6: (Color online) Same as Fig. but the solid curves (Fit II) are the best-fit results to 
reproduce the solar isotopic r-pattern for 69 < A < 124 with Tg = 1.5 and 10 20 < n n < 
10 22 ' 5 cm -3 , which are the typical conditions required by the N = 50 CWP nuclei. See text for 
details. 

To find the best match to the r-pattern in r-II stars, we consider a superposition of 
neutron densities in the two ranges adopted for Fits I and II. The results are shown as 
the solid curves in Fig. [71 It can be seen that fair agreement between the calculated and 
observed patterns is obtained for the FRDM, HFB-17, and WS* models. Note also that 
the calculated Te abundances (crosses) for all four models are consistent with the newly 
measured value for the metal-poor star BD +17°3248. However, the trough at Z — 66-70 
is clearly problematic for the RMF model. This may be caused by nuclear shape transition 
before the N = 126 closed neutron shell and the location of the transition region could have 
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been assigned incorrectly in the RMF model [50[. Discrepancies can also be seen for Ru, 
Rh, and Ag [Z = 44, 45, and 47, respectively) for all four models. This issue needs to be 
addressed by detailed considerations of the astrophysical environments for the r-process and 



alternative sources for the elements below the peak at A ~ 130 
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FIG. 7: (Color online) Best fits (solid curves) to the r-pattern in r-II stars using a superposition 
of neutron densities within the ranges 10 20 < n n < 10 22 5 cm -3 and 10 23 ' 5 < n n < 10 27 5 cnr 3 
for four nuclear mass models: (a) RMF, (b) HFB-17, (c) FRDM, and (d) WS*. The crosses give 
the calculated Te abundances, which are consistent with the measure value (filled square) for the 
metal-poor star BD +17°3248 (shifted according to the observed Te/Eu ratio). See text for details. 



B. Comparison with a non-solar-like r-pattern 

In contrast to the r-II stars, some metal-poor stars exhibit an r-pattern that is clearly 
different from the solar one. Prominent examples are the metal-poor stars HD 122563 65] 



and HD 88609 



661 ] . which have almost the same abundances for Cu (Z = 29) and heavier 



elements. The data on Sr [Z = 38) and heavier elements for HD 122563 are shown as the 
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filled circles in Fig. |HJ Relative to the solar r-pattern translated to pass through the Eu data 
(dashed curve), the elements below the peak at A ~ 130 in this star are grossly overabundant 
while those heavier than Eu are underabundant. It was argued that in addition to an r- 
process source for producing a solar-like r-pattern from the peak at A ~ 130 to that at 
A ~ 195, a very different source is required to explain the data for stars like HD 122563 
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FIG. 8: (Color online) Best fits (solid curves) to the non-solar- like r-pattern (filled circles) in the 
metal-poor star HD 122563 for four nuclear mass models: (a) RMF, (b) HFB-17, (c) FRDM, and 
(d) WS*. These results use a superposition of neutron densities within the range 

10 20.5 < 

n-n < 

10 25.0 

cm , which corresponds to the typical conditions required by the N = 82 CWP nuclei for 
Tg = 1.5. The dashed curve gives the solar r-pattern translated to pass through the filled circle for 
Eu (Z = 63). See text for details. 

Here we attempt to interpret the abundance pattern observed in HD 122563 using the 
classical r-process approach. We find that the conditions required by the N = 82 CWP 
nuclei can best reproduce this pattern while those required by the N = 50 and 126 CWP 
nuclei can not. The best-fit results use a superposition of ten neutron densities (equidistant 
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on log scale) within the range 10 20 ' 5 < n n < 10 25 cm -3 , which corresponds to the typical 
conditions required by the N = 82 CWP nuclei for Tg = 1.5 (see Fig. H]). These results are 
shown as the solid curves in Fig. [SJ It can be seen that an approximate overall match of 
the calculated with the observed abundances is obtained for all four nuclear mass models. 
Therefore, we suggest that it is plausible to account for the abundance pattern in stars like 
HD 122563 by an r-process operating under the conditions required by the N = 82 CWP 
nuclei. We note that Fig. [8] shows a clear difference in the calculated relative production of 
Te and Xe (Z = 52 and 54, respectively) between the WS* and the other three mass models: 
these two elements are produced in approximately equal amount for the WS* model but Xe 
is produced much more than Te for the other three models. Measurements of these two 
elements in HD 122563 would be extremely valuable in constraining nuclear mass models 
although they also represent a difficult challenge to spectroscopic observations. 



IV. DISCUSSION AND CONCLUSIONS 



We have explored the effects of four nuclear mass models (FRDM, WS*, HBF-17, and 
RMF) on the conditions required by r-process nucleosynthesis under the WP approximation. 
As discussed in Sec. [Ill the required T g -n n conditions are mostly determined by the two- 
neutron separation energies of the CWP nuclei with N = 50, 82, and 126 and of those 
nuclei around them. Figure [3] shows the dramatic effect of using the tabulated values in 
the latest atomic mass evaluation AME2011-preview when they are available to replace the 
masses predicted by models. As noted in Sec. Ill the tabulated masses of 76 Ni to 78 Ni 
and 78 Zn to 82 Zn play crucial roles in determining the conditions required by the N = 50 
CWP nuclei. However, the tabulated masses of 76 Ni to 78 Ni and 82 Zn are extrapolated 
rather than measured. To emphasize the effects of these masses on the conditions required 
by the N = 50 CWP nuclei, we first repeat the calculations of Sec. [TT] by varying the 



neutron separation energy of 78 Ni within the estimated uncertainty of 0.946 MeV 34| while 
keeping the other input the same as for Fig. [3b. The results are shown in Fig. [9K In 
comparison with Fig. [3b, the lower bound on the region of the required T 9 -n n conditions 
stays the same because this is determined by the two-neutron separation energy of 82 Zn 
(see Sec. Ill Bf ). which is not changed. Increasing the neutron separation energy of 78 Ni 
by 0.946 MeV raises the upper bound from the solid curve (upper bound in Fig. [3b) to 
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the dashed curve and decreasing this quantity by the same amount lowers it to the dotted 
curve. We then repeat the same calculations but vary the neutron separation energy of 82 Zn 
within the estimated uncertainty of 0.401 MeV {34J. The effects on the lower bound on 
the region of the required T 9 -n n conditions are shown in Fig. [§b. Note that if the neutron 
separation energy of 78 Ni were lower than its tabulated value by 0.946 MeV while that 
of 82 Zn were higher by 0.401 MeV, then it would be almost impossible to find any T$-n n 
conditions to accommodate all the N = 50 CWP nuclei. In any case, the significant effects 
of uncertainties in neutron separation energies on the required Tg-n n conditions shown in 
Fig. O clearly demonstrate the importance of precise mass measurements for 76 Ni to 78 Ni 
and 82 Zn. 




FIG. 9: (Color online) Effects of the uncertainty in the neutron separation energy S n for (a) 78 Ni 
and (b) 82 Zn on the Tg-n n conditions required by the N = 50 CWP nuclei. The solid curves are 
the same as those in Fig. [3b. The shaded regions show the effects on the required Tg-n n conditions 
when the corresponding S n values are varied within the estimated uncertainties. See text for 
details. 

As in the case of N = 50 CWP nuclei, we have also made a careful scan of the nuclear chart 
around the N = 82 CWP nuclei and explored the effects of those nuclei with experimentally 
unknown or poorly measured masses on the T$-n n conditions required by the N = 82 CWP 



21 



nuclei. When nuclear masses are not known experimentally, we have used the extrapolated 
masses and uncertainties as these have been proven to have a better predictive power than 
all available models 



67 



681 ] . We have identified 131 Cd and 132 Cd as additional key nuclei for 



precise mass measurements. 

In conclusion, we have estimated the temperature and neutron density conditions required 
for r-process nucleosynthesis under the WP approximation using four nuclear mass models 
augmented by the latest atomic mass evaluation AME2011-preview. We have shown that 
these conditions are mostly determined by the two-neutron separation energies of the N = 
50, 82, and 126 CWP nuclei and those around them. We have also identified some key nuclei 
including 76 Ni to 78 Ni, 82 Zn, 131 Cd, and 132 Cd for precise mass measurements at rare-isotope 
beam facilities. 

Based on the typical conditions required by the N = 50, 82, and 126 CWP nuclei shown 
in Fig. HJ we have performed r-process calculations in the classical approach to reproduce 
the r-pattern inferred for the solar system and those observed in metal-poor stars of the 
Milky Way halo. We have found that (1) the conditions required to produce the peak at 
A ~ 80 differ greatly from those required to produce the solar r-pattern from the peak 
at A ~ 130 to that at A ~ 195, which reflects that the T 9 -n n conditions required by the 
N = 50 and 126 CWP nuclei are very different; (2) the solar r-pattern from the peak at 
A ~ 130 to that at A ~ 195, which also closely describes the r-patterns in many metal-poor 
stars, can be reproduced under the conditions required by the N = 126 CWP nuclei, which 
has significant overlap with those required by the N = 82 CWP nuclei, thereby enabling 
coproduction of the peaks at A ~ 130 and 195; (3) it is plausible to explain the overall 
r-patterns in metal-poor r-II stars with a superposition of two sets of r-process conditions 
required by the N = 50 and 126 CWP nuclei, respectively; and (4) the non-solar-like r- 
pattern observed in metal-poor stars like HD 122563 can be accounted for by the r-process 
conditions required by the N = 82 CWP nuclei. 

We recognize that the classical r-process approach leaves out many important details, 
such as the time evolution of temperature and neutron density, the finite duration of the 
freeze-out, and the breakdown of (n, 7) (7, n) equilibrium during the freeze-out. However, 
it is very likely that the conditions derived here are valid approximation to those obtained 
immediately before the freeze-out in hot r-process environments where (n, 7) ^ (7, n) equi- 
librium can be achieved. We intend to carry out parametric studies of the r-process based 
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on more detailed and more realistic astrophysical models in the future, and will explore the 
effects of various nuclear input on such models. 
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